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I. INTRODUCTION 



The Kondo problem devoted to study the effect due to the exchange interaction between the 
impurity spin and the electron gas has played an important role in condensed matter physics since 
its discovery .The original treatments in Kondo problem the electron-electron interaction is 
discarded. This is reasonable in 3D where the interacting electron system can be described by 
Fermi liquid. Recently, much attention has been paid to the theory of the magnetic impurities in 
the Fermi liquid and Luttinger liquid, [Q, Q where the central scheme is the few impurity coupled 
with strongly-correlated electron system. Apart from the fundamental theoretical interests, it is 
remarkable that the physics implied here can be accessible experimentally. The recent advances in 
semiconductor technology enable to fabricate very narrow quantum wire which can be considered 
one-dimensional, and furnishes a real system of Luttinger liquid. Also edge states in a 2D electron 
gas for fractional quantum Hall effect can be considered as Luttinger liquid Q . Intense efforts and 
much progress has been made around the subjects from different approaches. Using bosonization 
and renormalization techniques, Kane and Fisher |^]studied transport of a ID interacting electron 
due to potential barriers. Their results triggered the study of the problem of local perturbations to 
Luttinger liquid and Kondo problem in Luttinger liquid. The Kondo problem in Luttinger liquid 
was considered by Lee and Toner |^ . They also performed the renormalization group calculation 
and found the crossover of the Kondo temperature from power law dependence on the Kondo cou- 
pling constant to an exponential one. Relying on poor man's scaling method, Frusaki and Nagaosa 
0] showed that the Kondo coupling flows to the strong-coupling regime not only for the antiferro- 
magnetic case but also for the ferromagnetic case. The boundary conformal field theory Q allows a 
classification of critical behavior for Luttinger liquid coupled to a magnetic impurity. It turns out 
that there are two possibilities, a local Fermi liquid with standard low-temperature thermodynam- 
ics or a non- Fermi liquid The non- Fermi liquid behavior is induced by the tunneling effect of 
conduction electrons through the impurity which depends only on the bulk properties but not on 
the details of the impurity B . Density matrix renormalization group calculation also supports the 
same conclusion . In addition the renormalization group flow diagram for parameters charac- 
terizing impurity is more complex and contains fixed points responsible for the low temperature 
behaviors when the potential of impurity is not strong |pT[| . 

Despite all important progress hitherto made, the problem of few impurities embedded in a 
strongly correlated ID electron system is still far from a complete understanding. We think that 
exact solutions of some integrable models on the subjects are useful from which one can expect to 
draw definite conclusions. Indeed Bediifig et al has thoroughly solved an integrable model with 
impurity coupled with t — J chain p^ . They introduce the impurity through a local vertices 
as in ||l^. The model introduced suffers the lack of backward scattering and the presence of 
redundant terms in the hamiltonian . Based on Kane and Fisher's observation we see it 
is advantageous to use open boundary problem with the impurities at open ends to study the 
problem of impurities coupled with strongly-correlated electron system. The program has been 
initiated for S~ interacting fermi system in for t — J model in and for Hubbard model in 
0]. 

The t — J model, is considered as one of the most fundamental models in strongly correlated 
electron system for its possible relevance for purely electronic mechanisms for high-Tc superconduc- 
tivity and heavy-fermion system. This model is obtained from the Hubbard model as an effective 
hamiltonian for the low-energy states in the strong- correlation limit. In this limit double occu- 
pancy of fermions is forbidden, leading to only three possible states at each lattice site for half 
spin. Currently, there is upsurge for its study. Very recently, the Luttinger liquid properties of the 
t — J model are discussed in Ref. By solving the functional relations, the finite-size corrections 
related io t — J model are dealt with for the open boundary conditions in Ref. jl^.The effects 
about an integrable impurity coupling to both spin and charge degrees of freedom are studied in a 
periodic t — J chain w hich we have mentioned above. The another generalization of the t — J 
model is given in Ref. ]19|] by using the one-parametric family of four-dimensional representations 
of f;^(2|l). It is also a kind of generalization of extended Hubbard model ]2C| ]. 

In this paper we expand the study of the Kondo problem in ID < — J model by exact solution 
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of open boundary Bethe ansatz. For this purpose we put two magnetic impurities in both sides of 
the open t ~ J model which is a typical situation for the one-dimensional systems with impurities. 
The coupling constants of the impurities with conduction electrons cover from negative infinity 
to positive infinity, which means that both the ferromagnetic Kondo effect and antiferromagnetic 
Kondo effect can be dealt with on the same setting. We then construct the hamiltonian for the 
system with magnetic impurities from t — J model. The integrability of this model ensures that 
both the Yang-Baxter equation and the reflecting Yang-Baxter equation are satisfied. By using 
the algebraic Bethe ansatz scheme for open boundary |21| we diagonalize the hamiltonian for the 
present system and obtain the Bethe ansatz equations. From which we derive the nonlinear integral 
equations governing the thermodynamic properties of the model for large system. The finite-size 
corrections for energy of ground-states in all cases can be calculated. 

The arrangement of the present article is as follows. In section 2 the constructed hamiltonian 
and its first quantization form are given explicitly. In section 3 the boundary matrix depending 
on the rapidity and spin of the particle is given and all possible integrable cases for the model 
are exhausted. The Bethe ansatz equations of the systems for all integrable cases are derived in 
section 4. The properties of the ground state for the cases other than that in [|l5|are discussed in 
great detail in section 5. In the final section the finite-size corrections of the ground-state energies 
for chosen cases are obtained. 



II. THE HAMILTONIAN OF THE MODEL 



Consider one-dimensional lattice with G sites , N electrons and two magnetic impurities at 
both ends. Due to a large on-site Coulomb repulsion there are at most one particle at one site. 
The dynamics of the system governed by a hamiltonian which we construct from the t — J model 
[22-27]. The conduction electrons can hop ( i ) between the neighbor sites. There are four types 
interactions in the model. A spin exchange interaction J and a charge interaction independent 
of the spin of strength V between the neighbor conduction electron; Kondo coupling Jq, Jband 
impurity potential interactions Vq, V^between the electron and impurities. The hamiltonian of the 
system has the form: 

G-1 G-1 G-1 

^ = E E (^/.c-.+i. + Q.) + J E • s,+i + ^ E ^^■"^■+1 

+ Ja^l ■ Sa + VaUi + JfeSc ' Sfe + VbTT-G, (1) 

where C^^lCja) is the creation (annihilation) operator of the conduction electron with spin a on 
the site j; Ja,b, Vafi are the Kondo coupling constants and the impurity potentials, respectively; 
Sj = i g., Cj'^acr,a'Cjcr' IS the spin operator of the conduction electron; nj = C^^Cj^ + C^^Cjiis 
the number operator of the conduction electron; G is the length ( or site number ) of the system. 
Some properties of the ground state for t=l, J=2, ^ = 3/2 have been reported in Ref. [T5[ |. 
Following Schultz's notation [^we write the translation operators T^: 

T^^{xi, ■■ ■ ,Xj, - ■■ ,xn) = '^{xi, ■ ■■ ,Xj ±1,- ■ • ,XAr), 

where • • • , Xj, • • • , xn) is the wave function of N conduction electrons. In first quantization 

and in appropriate energy units { t — 1) the hamiltonian of this system can be written down as 

N N 

^ = - E (^^ + + E + Kb,5.^, G + K,) (2) 

where the couplings are denoted by operators Kaj — Va — ^ + ^Paj and if^j = — ^ + ^Pbj 
with the permutation operators Pa(b),j between the spins of the conduction electron j and the 
impurities a, ( b). The operator Kj acts on the wave function as 
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JV 

Kj'^{xi,X2, ■ ■ ■ ,Xn) = ^Sx^.x, + lKij^{xi,X2, ■ ■ ■,xn), 
i=l 

where Kij = V — j + ^Pij describes the interactions between the conduction electrons with the 
permutation operator P^- permuting i -th and j— th electron in spin space. We will diagonalize 
the above hamiltonian in the following section. 



III. INTEGRABILITY CONDITIONS 



We write the wave function in region < xqi < xq2 < ■ ■ • < xqn < L — 1 as 



cri,cr2,---.<Tiv 



{XI,X2, ■ ■ ■ ,Xn) 



{rpQikpQi,rpQ2kpQ2, ■ ■ ■ ,rpQNkpQN) 

P ri,r2,---rH=±l 
N 

■ exp [i rpj k pj xj ] (3) 



where the coefhcients A 



crQi,aQ2,---,tTQN{rpQikpQi,rpQ2kpQ2, - ■ ■ ^rpQNkpQN) are also dependent 
on the spins of magnetic impurities which are suppressed for brevity,and sp — 1(— l),when the 
parity of P is even(odd), Sr — 11^=1 which r takes the value+lor —1. The boundary R matrix 
satisfies the reflecting Yang-Baxter equation: 



5'i2(A, Ai) R (A)S'i2(A, -m) R (m) (m)5'i2(A, -^i) R (A)S'i2(A, ^i), 

1 2 

where operators R (A) and R {fi) are defined as 

1 2 

R (A) = R{X) (g) idv2, R = ^dvi «> R{fi) 
for matrix R G End{V). S matrix satisfies the normal factorizablc condition: 
Suik, X)Si3{k, fi)S23{\, fi) = S'23(A,//)S'i3(A;,^)S'i2(fc, A). 



(4) 



(5) 



For convenience we set t = I. From reflecting Yang-Baxter equation and the form for S matrix, 
we refer that the boundary R matrix should have the form 



R — exp(i(^) 



q-iC-iP 
q + iC + iP' 



(6) 



where P is the permutation operator, q — ±icot|, ±itan| and C is the arbitrary constant. 
Putting Ka[b),j — m + IP , we have from eq. that 



1 k fc r 

(m - 1)^ - ^g^^2 _j^2l (g2 + 6-2 - 1) tan - + g \^{m + if - V 



^ 0. 



(7) 



This is the restriction imposed on coupling constants in order that our model (|l|) to be integrable. 
The details are as follows. 

In this case we know that the scattering matrix in the bulk can be written as: 

lPl2 



5l2(fcl, k2) 



2 cot 2 2 2 



i cot ^ - i cot ^ - i 



(8) 
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where P12 is the permutation operator between two electrons. The boundary R matrix at the left 
end of the chain takes the form: 

i?.(fc,,a,) = exp[^^„(fc,)]i^^li-:^^^i^. (9) 

The coupling constants Ja, Va at the left end of the chain are expressed in terms, of Ca 

J . ? V = ^^i^ (10) 

{2Ca T 1) i2Ca ± 3) ' '^'^ {2Ca T 1) (2a ± 3) ' ^''^ 



and 



Ja(cotf + 2za)exp(ifc,)+i[4+(4K- Ja)exp(iA:j)] ,,,, 
exp zo9a(/e.,) = ; . (11) 

Ja(cot ^ - 2iCa) exp(-i%) - z[4 + (4K - Ja) cxp(-^fc,)] 



The boundary R matrix at the right end of the chain has the form: 

1 cot ^ — iCh — iPi 

i?b(-fc,,(7,) = exp[-2^fc,(G + l)+^(^^,(fc,)]^ ^ ^. , , 

2 cot -J + iCfc + iPbj 

Similar relations exist for Jb, Vb and (pb{kj) by merely substituting indices a in ( [l^ ) and ( |ll| ) by 6. 



(12) 



A. J = -2 V = i 



The boundary R matrices have the forms: 



^ tan -1^ 4" ^Cfl -t- ^-^aj 



Ra{kj,a,) = eM^Mkj)]i ^ — ^, (13) 

2 tan iCa *^aj 



i tan — + iCi + iP;, 

i?b(-fc„a,) = eM-2ik,iG + l) + iipb{kj^ ^ . ' (14) 

i tan ^ — lOf, — ziijj 

ipa{kj) and ipb{kj) are the same as in the proceeding. Now the coupling constants should be written 
in terms of the arbitrary parameter Co in the form 

J _ 8 y _ 4Cg - 3 

" (2Ca T 1) (2Ca ± 3) ' " (2Ca T 1) (2Ca ±3) ^ 

Jh. Vb have the same expressions except with the substitution of indices a by b. Correspondingly, 
the scattering matrix S in the bulk for two conduction electrons is 



i tan 4f - i tan ^ + iPj 



12 



2 tan "2 tan -h z 
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B. J = 2 V = I 



In this case the dependence of couphng constants on parameter Catakes the form 

j= § " ^ (17) 

Jb. Vb have the same expressions by the substituting of indices a by b. The scattering matrix in 
the bulk is 

i tan ^ - i tan ^ - iP^ , , 

5i2(fci,fc2) = \^ \^ k.J - (18) 

^ tan ~2~ ^2 tan ~2~ ~t i 



The boundary R matrices are 



^ tan — iCa *-fai 



Ra{k,,aj) = eM^Mkj)]Y- ^ — ^, (19) 



- tan — — iCi — iP?, ■ 

Rb{-k^,a,) ^ cxp[^2ik^{G + I) + iMkj)]j ^ ^ -■ (20) 

2 tan + iCb + iPbj 



C. J = -2 V = 
The couphng constants have the forms 

(2CaTl)(2a±3)' " (2aTl)(2Ca±3)" ^ ' 

Jb, Vb have the same expressions by the substituting of indices a by b. The scattering matrix in 
the bulk is 

W/^ A-)- jcotf -Icotf +»Pi2 
The boundary R matrices are 

~r cot "TT" H~ iC'n ~\~ 1 1^n o 

Raik^a,) = exp[z(p,(fc,)] ^ ,^_-/„-p ' ^^3) 

2 cot 2 ^ o-j 

1 cot — + iCh + ■ 

i?b(-fc„(7,) = exph2*fc,(G + 1) + *^f,(%)]^ ^ ^, (24) 

2 cot -J — iCfc — iPbj 

The expressions for boundary matrices depending on both the moment of the particle and the spin 
of the electron are new. The expressions of S matrix in the bulk have been obtained before in Ref. 
p3],but they are different from ours. 
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IV. BETHE ANSATZ EQUATIONS 



By using the standard Bethe ansatz procedure, we can diagonalize the hamihonian (|l|) ||2l[| and 
obtain the following Bethe Ansatz equations. When J = 2 and V — setting 

i cot ^ - i cot % - iP,o , ^ 

^^o(k,,ko) = \ \ (25) 

^ cot - ^ cot ^ - I 

(, . 1 cot ^ - ^ cot ^ - tPjN+1 

W) - ' (26) 

where cot ^ — 2iCq, cot '^"^^ — 2iCb, Pjo = Paj,PjN+i = Pbj, we can write down the boundary 
R matrices as the forms: 

R.{k„a,) = exp [^Mk,)] j | ' ' ' ^^^'T' ^ , , (27) 

i cot ^ + iCa + I Sja{-kj,ko) 

1 k ■ 



i cot ^ + iCb + I SjN+i[~kj,kN+i) 

Define 

r(A) - Srj{X)SroWSrl(X) ■ ■ ■ Srj-l{X)Srj + l{X) ■ ■ ■ SrN{X)SrN+l{X) (29) 



with 



A - i cot |- - iP, 
A - i cot ^ - i 



SriiX)= ^ \ \ ; = 0,l,---,iV+l. (30) 

i cot f 



We get the equation 

r.[T(A)r-(-A)]l ^, ^ 



2i — cot ^ \ cot ^ + iCa + * 5 cot ^ + iCh + i 
i — cot ^ cot — zCa — i 5 cot — iCb — i 

•exp [-i(fa{kj) - ifb{kj) + 2ikj[G + 1)] $ (31) 
where $ is the eigenstate of the system. Then the Bethe ansatz equations can be expressed as 

exp [2ik,{G + 1) - i^a{kj) - iipb{kj)] \ 1 :| 1 : 

^ cot — iCa — I cot -J — iCb — I 



T-r ^ cot ^ - A/3 + ^ i cot ^ + A/3 + ^ , , , ^ 

p=l 5 cot ^ - A;3 - H cot ^ + A;3 - t 



(A„ + §)2 + C2 (A„ + §)2 + C2 -jX- A„ - i cot -I + I Aa + i COt ^ 



n 



(A„ - 1)2 + C2 (A, - 1)2 + C2 f \ A„ - i cot I - i A, + i cot I 
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A'l 

- n («-l-2,.,M). (33) 

M is the number of down spins and N is the number of the electrons. The function ip is denoted 
by expression ( ^l]). Similarly, when J = —2 and — ^, we can write down the Bethe Ansatz 
equations as the forms: 



exp [2ikj{G + 1) - i^a{kj) - iLpb{kj)] - 



2 tan ^ ^ tan ^Cf^ 



tan + iCa + 2 i tan ^ + zCb + i 



i tan -rf — Afl — i tan -rf + Aa — -I , , , , 

2 2 ^ 22 2 ^ 2^ (j = 1,2, •••,7V) (34) 

;3=i ^tanf -A/3 + I itanf + A/3 + f 



(A„+|)2+C2 (A„ + i)2 + c2 .^ A„ -itan| + i A„ + itanf + 



(Aa 2)'^ ^" C'^ (Aq 2)'^ + C'ft ;_-|^ Aq ^ tan ^ Aq + ^ tan ^ — 



M 



Xa — Xfj + i Xa + + i 



n :"'7" c-Tr" (a=i,2,^^^,M). (35) 

\a ~ Xjj — I \a + Xp — I 



When J = 2 and F = |, we have that 

s , , , tan + iCa + j ^ tan + iCb 
exp {2U,,{G + 1) - z^„(fc,) - z^b(fc,)] f— ^ . ° . ^ ^ 

2" tS^n ~2~ ^2^ Lclll ~2~ 



^ i tan ^ - i tan ^ + z i tan ^ + i tan ^ + I 



n T 

i=i{i^j) 2 ^^"^ ^ ~ 1 ^^"^ "2" ~ * § ^a'^ "2" + 1 "2 



7^ tan ^ - An - i i tan ^ + A^ 



(Act + 5)^ + (Aq + ^)^ + -j-j- Aq — ^ tan + ^ A^ + ^ tan ^ 



(Aq, — ^)^ + (Aq. — ^)^ + Aq — ^ tan — | Aq + ^ tan ^ ~ 

= n ""'""^'""t""^' (a=l,2,^^^,M). 

A„ — A3 — z A„ + A3 — z 



When J = — 2 and F = -|, we get that 



i cot ^ — iCa — j ^ cot ^ — iCf, 



exp [2tk,{G + 1) - z¥Pa(fc,) - i^k/)] I , fc, , . „ , . 1 , fc, , , . 

i cot + iCa + I 5 cot ^ + iCb + I 

-Q i cot ^ - i cot ^ - i i cot ^ + i cot ^ - i 



(37) 
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1 cot ^ - - I i cot ^ + 



n 



/3=1 5 cot ^ - A/3 + I i cot ^ + A;; 



(j = 1,2, •••,7V), (38) 



(Ag + |)^ + C,^ (A„ + |)^ + -jL A^ - ^ cot f + f A„ + 1 cot f + f 



(A„ - f )2 + CI {K - 1)2 + Cl A, - i cot f - i A„ + i cot 



2 



M 

Aa — A/3 + i Aq, + A/3 + i 



n 



Arv — A/3 — Z A^ + A/3 



(a=l,2,.^^,M). (39) 



Here the function ipa{kjj is expressed by equation (^jj) and ipi,{kj) has the same expression as 
relation (|l^) with the substitution of index a by b. M is the number of down spins and N is the 
number of the electrons. It should be noted that in the above Bethe Ansatz equations we have 
choose the boundary R matrices as the form as the ones in Section 3, where the boundary matrices 
depend on the spin parameter. If the R matrix is dependent on the spin of the electron only at 
one end of the chain, for example, we denote by Ri,{kj,aj) the boundary matrix at right end of 
the chain independent on the spin aj . The Bethe Ansatz equations for J = 2, V = — ^ take the 
form: 

exp [-iipajkj)] 2 cot ^ + iCg + I ^ ^ cot ^ - A/3 + § g cot ^ + A/i + § 

Rb{~kj,dj) 1 cot ^-zCa-z jJiicot^-A/3-f icot^ + A/3-|' 

(j = 1,2,. ••,]¥) 
(Aq + ^)2 + -j-j- Aq — ^ cot + ^ Aq + ^ cot + f 



(A„ - 1)2 + C2 f^-^ A„ ~ icot^ - f A„ + icot^ 



2 



M 

n 

/3=1(/3#q) 



Aa — A/3 + z Aq + A/3 + i 
Aq — A/3 — 2 A„ + A/3 — i ' 



(a = l,2,.^^,M) (41) 



Similarly, when the boundary matrix at the left end of the spin is independent on the spin of the 
electron, denoted by Ra{kj we have that 

exp [2ikj{G + 1) - i(pb{kj)] ^ cot ^ + iCt + i 
Ra{kj,dj) ^cot^ -iC'b-i 



^ 1 cot ^ - A/3 + ^ i cot ^ + A/3 + 



n 



2 2 'P ' 2 2 ^" 2 ' "P ' 2 

1 cot ^ - Afl - H cot ^ + A,3 - ^ 



{j^l,2,---,N) (42) 



13=1 2 2 '~P 2 2 2 



(Aa + 5)2 + -i r Aq — ^ cot + ^ Aq + ^ cot ^ 
(Aq - f )2 + C2 11 AQ-icot^-§AQ + icot^ 



Af . . . . . 

n ^''^^^^t'^'^^ («=i,2,.^^,m) (43) 

Aq — A/3 — z Aq + A/3 — I 

/3=1(/35^q) " P a p 
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where the number of down spins should less than N + 2 and N is the number of the conduction 
electrons in the system. Furthermore, we get that 

exp [-i'Paikj)] ^ tan ^ - iCq - i ^ jl ^ tan ^ - Ag - | ^ tan ^ + Ag - | 
Rb{-kj,aj) ^t&n!f +iCa+i j-i^ i tan ^ - A^ + f i tan ^ + A;3 + f ' 

{j = l,2,...,N) 

(Aa + + C'a 'YJ^ '^'^ ^ \ tBX^ ^ + 5 •^a + \ tan ^ + 5 



(Aa - 1)2 + C2 ^^-L A„ - itanf - i A„ + itanf - 



M 



n :""'"^'^t^'^' (a = l,2,...,M) (45) 

Act — Afl — z Aq, + A/3 — I 



and 



exp [2ikj{G + 1) - iiph{kj)] \ tan ^ - iCb - i 
Ro.{kj,dj) 1 tan ^ + iCb + i 



i tan ^ - A/3 — ^ i tan ^ + A/3 — i , ^ , ^ 

n L„„/:^ ^ , , (i = l-2,---,A^) (46) 



ii5tan^-A/3 + titan^+A^ + i 



(A„ + ■^)2 + -j-j- Ad — i tan + ^ A„ + i tan ^■ 



{}^a 2)'^ ^6 ;=i '^a ~ 5 ^an ^ ~ Aq- + ^ tan 



Aa - A/3 + Z Aa + A/3 + i 



n 7T7 (a = l,2,..-,M) (47) 

Aa - A/3 - I Ac + A/3 - Z 



for the case of J = —2, V = ^. 

exp [— iVa(fcj)] 5 tan ^ + + i YJ \ ^ ~ | tan \ +i\ tan ^ + | tan ^ + i 
Rb{-kj,dj) 1 tan ^ - iC„ - i i tan ^ - i tan f - i i tan ^ + i tan f - i 

^^ jt.„|-V + Hta„| + A, + j ^ (, = l,2,..,iV,, ,48, 
/3=i 2 ^ - A/3 - 2 2 ^ + A^ - 2 



(Aa + ^)^ + C"^ -I r Aa — \ tan ^ + ^ Aa + ^ tan 
(Aa - 1)2 + C2 11 Aa - i tan f - i Aa + i tan 



n ^"'" + '^tx'^' (a = l,2,...,M) (49) 

Aa — A/3 — Z Aa + A/3 — I 
/3=l(/3^a) " a I p 
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and 



1 k ■ 

exp [2ikj{G + 1) — iipb{kj)] ^ tan ^ + iC'b + i 
Ra{kj,dj) ^tan^-iCb-i 

^ ^ tan ^ tan ~t~ ^ ^ tan -t- ^ tan -t- z 

2 t^"^ '2 ~ 2 "2" ~ * 2 t^"^ "2" + 2 ^^"^ "2" ~ * 

n''^ i tan ^ — Af) + ^ i tan + + -I . s . ^ 

L 4 X (i = 1.2,---,A^), (50) 

p=l 2 T ~ "^/? ~ 2 2 "2" "^/? ~ 2 

('^a ~l~ ■^)^ ~^ ^6 j^j" "^ci: ^ tan 4~ ^ H- ^ tan ~l~ 



(Aa - 1)2 + C2 f \ A„ - i tan f - i Aa + i tan 



1 kl 



M 



Aa - A/3 + i Aa + A/3 + i 

/3=l(/3#a) 



= n r~r 717 (a=i,2,...,M) (5i) 

Aq,- A/3 -I Ac + A/3- z 



for the case oi J = 2, V = | when boundary matrix only at one end of the chain rely on the 
spin parameter of the electron. Finally, for the case of J = — 2 and V = the Bethe Ansatz 
equations take the form: 

exp [-i(pa{kj)] I cot ^ - iCa — i ^ cot ^ - ^ cot ^ - i i cot ^ + i cot ^ - i 

Rb{-kj,dj) i cot ^ + iCa + i i cot ^ - i cot f + i i cot ^ + i cot f + i 

n''^ i cot %- - Afl — ^ i cot + A/3 - ^ , , , , 

J, ' ' (J = l,2,---,n (52) 

/3=1 2 cot ^ - A/3 + I 2 cot f + A/3 + f 
(Aa ~t~ ^)2 ~t~ C'q j— j Aq; ^ cot "2" ^ Aq + ^ cot "2" ^" ^ 



(Aa 5)2 + C'a Aa — ^ COt ^ ~ 5 Aq, + ^ COt ^ ~ 5 



M 



Aa - A/3 + i Aa + A/3 + i 

/3=l(/3#c<) 



= n r"7".r:7". (a=i,2,.--,M) (53) 

■^-»- Aa - A/3 - Z Aa + A/3 - I 



and 



exp [likj{G + 1) - iipbikj)] ^ cot ^ - iCb - i 
Ra{kj,dj) icot^+iCh + i 

-Q i cot ^ - i cot ^ - i i cot ^ + i cot ^ - i 
i=i((^j) 5 cot -f — \ cot ^ + « 5 cot ^ + 5 cot ^ + i 



^J^lcot^-A/3-Hcot^ + A/3-i (. = 12... TV) (54) 
ii icot^-A/3 + iicot^+A/3 + i' 
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{^a H" H~ J^J" — ^ cot "2" ^ '^ct ^ C*-*^ ~2 



(Aa - + C2 ^^-^ A„ _ i cot f - f A„ + i cot f 



^ X -X,+^X +X, + ^ ^ 3, • • ■ , M). (55) 

Ra{kj,aj) and Rb{—kj,aj) denote that the boundary matrices at the left and the right ends of 
the system are independent on the spin aj respectively. Notice that the number of down spins is 
less than iV + 2 for the system with N conduction electrons. In the following section, we focus the 
discussions on the system with the boundary matrices depending on the spins of the electrons at 
both ends of the chain. Set 

1 (4C2 - 3) cos fc - 4Cg + 5 + 4iCa sin k 
' ^ t ^ (4C2 - 3) cos A: - 4C2 + 5 - AiCa sinfc ' 

1 (4C2-3) cos fc-4C2 + 5 + 4^a sinfc 
' " I (4C2 - 3) cos fc - 4C2 + 5 - 4iCb sin fc ' 



From relation (|ll|) and 

. ^, Jbicot!^ + 2iCb)exp{ikj) + i['i+{Wb- Jb)cxpiikj)] 

exph(pb[k.j)\ — 7-^ (57) 

Jb(cot ^ - 2iCb) exp(-ifcj) - ?;[4 + {Wb - Jb) exp(-ifcj-)] 

we get the following expressions. When J = 2 and V = we have that 

8 T/ _ 3-4C2 



exp [Z^.(fc)] = { '"P^f ) f - (2C.-l)(2C„+3) . - (2C.-l)(2C„+3) . (53) 

CXp[lfc + lOa(k)\ lOr - ^(2Ca + l)(2C„~3)' ~ (2Ca + l)(2C„-3) 



When J = —2 and y = ^, we have that 



2 ' 

4C^ 



exp(ifc) for Ja — /2r7„-lV2(7„+31 ' 



exp [^<^,(fc)] = _ . _ '^'^ - (2C„-l)(2C„+3) - - - (2C„-l)j2_C.+3) (59) 

[ exp[zfc i6'a(7r fc)]for Ja — (2C„ + l)(2C„-3) ' ~ (2C„+l)'(2Ca-3) 



When J = 2 and F = |, we have that 



„^„r,-,. _ / exp[2fc + z6'a(7r fc)]for Ja - (2C„-l)(2C„+3) ' " (2Ca-l}(2C<.+3) 

I -exp(lfc) for Ja = (2C„4-1U2C„-31 '^° ^ 



(60) 



(2Ca + l)(2Ca-3) ' '^'^ (2Ca + l)(2C„-3) 



When J — —2 and y = — |, we have that 



exp [Z^aik)] = I '^''Pf,'^ '^''^ " (2C.-l)(2C„+3) ' - (2C„-^lK2C„+3) ^ (g^) 

[ exp(ifc) for Jq = (2Ca + l)(2Ca-3) ' (2Ca + l)(2C„-3) 

The expressions of exp[i(^h(fc)] can be obtained by substituting the index a of 6 in the above 
relations. Then, without loss any generalization, we can choose that 

J- ? V 



(2a-l)(2Ca + 3)' (2Ca - 1) (2a + 3) ' 
J = 8 y ^ 3-4C2 

' (2^-1) (2a + 3)' " (2^-1) (2a + 3) ^ ' 
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for J = 2 and V = The Bethe Ansatz equations take the forms as 



n 



(A, - i)2 + d (A„ - i)2 + Cl f^^ Xo^-qi-^Xo^+qi-^ 

Aq — A/5 — z Aq + Afl — I 

where = i cot For the case of J = — 2 and V = ^, the Bethe Ansatz equations have also the 
forms ( |63| ) and (|6^ ) with qj = tan and the Kondo couphng constants should be 

4C2 - 3 



°" (2a-l)(2Ca + 3)' (2Ca - 1) (2a + 3) ' 
8 _ 4CI-3 

" (2a-l)(2Cb + 3)' " (2a-l)(2Cfc + 3)- ^^^^ 



If J = 2 and V = ^, we choose that 



4C2 - 7 



(2Ca + 1) (2Ca - 3) ' (2^ + 1) (2a -3)' 
8 4:C^ — 7 

" " (2a + 1) {2Cb - 3) ' " (2Cb + 1) {2Cb - 3) ' ^^^^ 



and the Bethe Ansatz equations are 



N 



exp(2zfc,G)^^±^(^^^±^(^ n (67) 



9j-»(C„+l) 9j-^(C,+l) -qi-iqj+qi-i 



-Q - A^ + 2 + A^ + 



,1 1j - - f 9i + 

1 ^ . They are alsc 

V = — I with qj = — i cot ^ and 



and relation ( |6^ ) with = ^tan-^. They are also the Bethe Ansatz equations for J — —2 and 



(2a + 1) (2a -3)' (2a + 1) (2a -3)' 
' (2a + 1) (2a -3)' " (2a + 1) (2a -3)- ^ ' 
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V. GROUND STATE 



In this paper we restrict the discussions of the properties of ground state to the case of J = ±2 
andF = TiThe case of J = 2 and V = ^ were studied in . The eigenvalue of the hamiltonian 
is 

E = T2N±J2-^ (69) 

J = l 'i j ^ 4 

for J = 2, F = — ^ with qj = i cot ^ and J — —2, V = ^ with = — i tan respectively. They 
satisfy the Bethe Ansatz equations (^3|) and ( |6^ ) from which the integral equations are derived. 

A. Integral Equations 

Following we introduce the notation 

, ■, X + i 
e{x) = :. 

Then, from relations (|6^) and (^^we get 

e I e (^^^] e(2g,)2G ^ J| g (2g, - 2X0) e {2q, + 2\p) , (70) 

(3=1 



l + Ca Vl + Cb 



N M 

■We{2K-2qi)e{2K + 2qi)= e (Aa - A^) e (Aa + A^;) , (71) 

where j — 1, 2, • • • , A^; a = 1, 2, • • • , Af and e(±oo) = 1. Considering that the parameter qj can take 
complex values, the general structure for {qj}j=i^2,...,N should be consisting of Af pairs of q^ — 
A„ ± I + 0(exp(-(5G),a = I,- - •,A//' and M" pairs of g = -A„ ± | + 0(exp(-5G), A„ G {A^} 
and remaining N — 2{M' + M") non-pairing q'^ s. To be more precise, we use 

Q = [q, I J = l,2,...,iV} - A'UA"ur, 

where 

X' = {q^ = K±^- + 0(exp(-^G) I a = 1, ■ • • , A/'}, 

X" = {qt = -K±'^ + 0(exp(-(5G) I A„ e {A^}, a = 1, • • ■ , Af"}, (72) 

r = g\(A'uA"). 

Obviously, the non-pairing satisfies equation ([70|) with j = l,2,---,iV — 2(Af' + Af"). When 



qj e X', from equation (70), we have 



Aq. A / Aq \ / Aq, \ / A, 



2G 



14 



M 

= e{2q+ ^ 2Xa)e{2q- ~ 2Xa)e{2Xa) [] e(A„ - A/3)e(A„ + A^), a = 1, 2, • • ■ , M'. (73) 

/3=l,(/3#a) 

From equation (j?^) we have 

T^) e (^) e (^) e e(2A.) 

,2 — CaJ \i2+Ca J — CbJ + Cb J 

■ W e{2K-2qi)e{2X^ + 2qi) \{ e(A^ - Aa)e(-A^ - A^,) 
1=1 p=i 

= ei2q+ - 2Xa)ei2q- - 2Xa) (74) 

where M_ = M'+M" and A^ ( /3 = 1, 2, • • • , M—MS) are the parameters describing the down spins 
but having no contributions to the bound states. With the help of the above relation, equation 
(73) becomes 

e ( e ( e ( e (A.)^« 

Ar-2M_ M' 

^e{2Xcf n e(2Aa - 2g0e(2A„ + 2g0 \[ e(Aa - A^)e(A„ + A/3) 

1=1 /9=l(/3#a) 

M" 

■\{e{X^-XfMK+Xf3), a = l,---,M'. (75) 

13=1 

Similarly, when qj G X",we get the following equation 

A \ /a \ /a \ /a \/~\2G 



\Ca + i r[ca-u\Cb + i r\Cb 



= e (2A„) W e(2A„ - 2gj)e(2A„ + 2qi) [| e(A„ - Xp)e{X^ + Xp) 

1=1 /3=1 

M" 

• Y[ e{X„-Xp)eiX^ + Xp), a=l,---,M" (76) 

/3=1(/3#q) 

The two equations (eell) and (eel2) can be combined in a single equation. 

Xa \ I Aa \ / Aq \ / Aq 



N-2M 

^2 



r ) e(Aa)^'^ 

5 / 



(2A„)' [| e(2A„ - 2g0e(2A„ + 2g0 



e 1^ 

M_ 

■ W e(A„-A^)e(A„+A/^i), a = l,2,...,M_, (77) 

/3=1(/3#q) 
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with the new Aq defined by 

Aq = 



Aq when a = 1, 2, Af ' 

Am'-o when a ^ M' + 1, M' + 2, M_ 



The parameters Aq (a = 1, 2, • ■ • , Af — in view of ([74[), satisfy 



2/ '2/ 1=1 



Aq \ / Aq 



Ca 2/ V'-^'' 2/ /3=1(/3^Q) 



The non-pairing (i,e. qj G y)satisfies 



^ {c^rd ' (c^) ^ ri ' (29, - 2A^) e {2q, + 2A^) 

M-Af_ 

e (2qj - 2A^) e (2^^ + 2A 



where j = 1, 2, • • ■ , iV - 2M_ and e(±oo) = 1. Setting 

e{x) = 2 tan~^ X, -tt < 9 < n, 

we have 

e{x) = exp[i(7r — 9{x))]. 
The logarithms of the equations (|7^), (|7|) and give, respectively, 

N-2M- 

= 4^Jq+0(2Aq)+ Y1 [^(2AQ-2gO+^(2AQ + 2g,)] 
+ ^[0(Aq-A0) + 0(Aq + A^)] 

/3=1 

with a = 1, 2, • ■ • , M_ and integers or half-integer Jq; 

7V-2A/_ 

M2\a - 2qi) + e{2K + 2(7;)] 



Ca + \] \ Cfc + i 



1=1 



47rjQ - 6I(2Aq) + 9 



-^Q \ n I 



CI I I 1 

a ~ 2 J \ b ~ 2 



M-A//_ 

■ 51 [^(Aa-A/3) + 0(AQ + A^)] 

/3=1 
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with a = 1,2, ■ ■ ■ , M — A/_ and integers or half- integer Jq,; 



Af_ M-M_ 

= 47r/, +^[0(2gj-2A/5) + 0(2g, +2A;3)]+ ^ [0(2?^ - 2A;3) + e(2g, + 2A^)] 

13=1 13=1 

with j = 1, 2, • ■ • , iV — 2i\f_and integers or half-integer Ij. By setting 



-^9 [k {x + c)] = 2TTa + ) ' 



the equations (18(1), (J81|) and (82) can be changed into the forms 



+ 2Ga (A„, 



dJa 
dXa 



a A„ - g;, - + a A„ -I- gf; 



- ^ [a (Aa - \f3, 1) + a (Aa + A/3, 1)] 
/3=1 



with a = 1,2, • • • ,M_; 



a ( Aa, Ca + - j + a f Ac, Cb + - 



E 



a ( Aq - g;, - j + a ( Aa + g;, - 



M-M- 



E [a pa 



/3=1 



A;3, 1 1 + a I A„ + Afl, 1 



V + -A/3, 



with a = 1,2,- • •,M - M_; 



a (gj, Ca + l) + a {q^,Cb + 1) + 2Ga ( qj, - 



13=1 



a ( 9j - A/3,i j + a Uj + A/3, ^ 



E 

/3=1 



a ( 9^ - A/3,^ j +a Uj + A/3,^ 



with j = 1, 2, • • • , iV - 2M_ . We define that 
J(A) ^ OiX) + 
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N-2M- M_ 



1=1 13=1 



N-2M- M-M- 



+ ^\ E mX-2qi)+0{2X + 2qi)]- ^ ie{X -X^) + 9(X + Xp)]\ , 



1=1 13=1 



h{q)^e{2q) + ^[0(-^]+i 



2G i \Ca + lJ \Cb + l 



M_ M-M- 



i ^ [0(2g - 2Xp) + ^(2g + 2A;3)] + ^ [0(2g - 2A^) + ^(2g + 2X0)] \ . 



2^ 

Then, the holes of A, A and q are defined as the sohitions of 

Gj{X) = 27r X (omitted J) , 



Gj(x) = 27r X (^omitted j) 



G/?,((7) = 27r X (omitted /) . 

By taking the thermodynamic hmits, we introduce the distribution functions 

X^a{\) \ (a\X) 
q ^ Q(q) hole. _ 

A -5(A) J {a\X) 



So we have that 



^=27r(a(A)+a^A)), 

'^J^ = 2n{p{q) + e\q)), 

^=2.(5(A)+a^A)). 
Therefore, the integral equations can be written down as 



2a(A,l) + i 



a[X,Ca + \] +aix,Ca-]:] + a( A,C6 + ^ ) +a(x,Cb-]: 



2 2 \ 2 \ 2 
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a [\, + 2cr(A) + 2(T^{X) + j d\'a{\')[a{\ - A', 1) + a(A + A', 1) 
dqp{q) 



a ( ^ - <?, ^] + a ^A + g, i 



2a [q. ^ + ^[a{q, + 1) + a{q, C^, + 1)] = 2p{q) + 2/(9) 



(92) 



d\a(\) 
d\d(X) 



a ( 9 - A, ^ j + a + A, i 
a ( 9- A, ^ j + a U + A, ^ 



(93) 



1 

G 



a(A,- 



a I A, - - 



a ( A - i j + a ^A + g, i 



2ct(A) + 2ct''(A)+ / dA'CT(A')[a(A- A',1) +a(A + A',1)], 



(94) 



where a(A, 77) = ?7/[7r(A^ + 77^)] with the arbitrary parameter 77. The terms with factors l/G in the 
above three equations describe the finite-size corrections of the system. 



B. Properties of Ground State 

For the system with N electrons, by using the distributed functions ct(A), 17(A) and p{q), the 
particle number and magnetization per unit length are given by 

N 
G ' 



dqp{q) + 2 / dAcr(A) 



dqp{q) — I dX5[\) 



G 2 

The energies per unit length have the forms as 

E IN f / 1 A f 

- = -— + 27: j dqp{q)a iq, ^]+'^'^ j d\a{\)a{\ 1) 

for the case ofJ = 2, y = — 1/2 and 

§ = ^ - 2^ y dqp{q)a (^q, ^ - 2vr ^ dAa(A)a(A, 1) 

for the case ofJ = -2,y = l/2. The relations ^ and (H) become as 

2a(A, 1) = 2a{\) + 2cr''(A) + / d\'a{\')[a{\ - A', 1) + a(A + A', 1) 



(95) 



(96) 



(97) 
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« ( ^ - 9, + a + i 



2a q, 



f J dqp{q) 
2p{q) + 2p^{q)+ I dXa{X) 



(98) 



a ( <? - A, ^ j + a U + A, i 



dXd{X) 



a ( 9- A, ^ j +a U + A, i 



(99) 



dqp{q) 



a[X-q,^j +a(x + q,^ 



if we set G 



= 2d{X) + 2d''{X) + J dX'd{y)[a{X - A', 1) + a(A + A', 1)], 
-oo. By Fourier transformation of equation ( |9^) we get that 
dqp{q) + 2 [ dXa{X) + I dX<7'\X) ^ 1, 



(100) 



which gives that N/G = 1 — J dXa^{X). Owing to <J^{X) > 0, we have that N < G, which coincides 
with the single-occupancy of every site. We assume that there is one particle per lattice site, that 
is, N/G = 1. Then we have cr'^{X) = 0. Now we consider the case of nonmagnetic p{q) = 0. The 
relation (|9q) turns into 



a(A, 1) = ct(A) + / dX'a{X')a{X - A', 1). 



(101) 



By Fourier transformation of equation ( 100 ) we have that a{X) — f?''(A) = 0, which means that 
Sz/G = and the system is nonmagnetic . From the above relation, we have that 



a{X) 



1 

2^ 



2 cosh I 



-div. 



(102) 



The interesting thing is that the above expression is exactly same as the integrable narrow-band 
model with periodic boundary obtained by Schlottmann p3] . In this way, relation (B^) reduces to 



and it gives that 



(q, = P'\q) + I dXa{X)a (^q ^ A, 1^ , 

i sec h \TTq\ , for q ^ 



271 



+ 00 



oo 2 cosh f 



diij = 



i, for q 



(103) 



The number M of the down spins is equal to G/2. The ground-state energy is E/G — —2 In 2 for 
J — 2, V = —1/2, which has the same value as the one in the periodic boundary condition It 
is due to that the impurities located at the both ends cause only the finite-size correction of the 
ground state energy. For the case of J = —2, V = 1/2, corresponding to the ferrimagnetic state, 
we have that 

1 k 



Piq) = - 



TT q- 



1 ' 

4 



a^A) = 



1 1 



TT A2 + r 



a{X) = 0, 



(104) 



by taking into account of cr(A) = o'''(A) = 0. Then we have E/G = 0. 



20 



VI. FINITE-SIZE CORRECTION OF THE GROUND STATE 



We assume that the distribution functions (t(A), p{q) and ^{X) are even functions about param- 
eters A, q and A, respectively. Then we have the following equations: 



a(A,l) + ^ 



a{ X,Ca + ^j + a(x, Ca - + a (x,Cb + + a (x, Ct - ^ 



a (^X, + a(A) + a'^(A) + J dX'a{X')a{X - A', 1) + J dqp{q)a (^X - 9, ^) , (105) 



1 /, 1 
2G 



a(g,^j +^Hq,Ca + l) + a{q,Cb + l) 



P{q) + P'('Z) + J dXa{X)a (^q-X,^^+ J dXd(X)a (^q - A, , 



(106) 



1 

2G 



a ( A, i j +a(x,Ca + + a (x,Cb + - a (x,Ca -^j -a(x,Cb-^ 



+ / dqp{q)a [X-q, 



1 



^d(X)+d''(X)+ / dA'CT(A')a(A- A',1), 



(107) 



from equations (|9^), (^) and (|94|), where a(A, 77) = r]/[n{X'^+r]'^)] with the arbitrary real parameter 
rj. The terms with factors l/G in the above three equations describe the finite-size corrections of 
the system. The energies of the system can be described by 



E 2N 



dqp{q)a q 



1 



dACT(A)a(A, 1) 



for J — ±2, V — t1/2, respectively. Setting 

Sri = sign{n) = 

we have that 



1, r;>0 

-1, ^<o, 

0, r; = 



a(w,r/) = S'^exp(- \ujri\). 



By Fourier transformation of equation (105), we have 



1 

2G 



for N/G = 1. By letting 



6(A)=a( A,a-f + aUCa-ij + ^ 1 +a\X,Cb ^ 



(108) 



(109) 



(110) 



a X 



(111) 
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we have that 



dqp{q)a (^q, ^) + / dXa{X)a (A, 1) - a(0, 1) + ^ - ^(0) - 'j\0). 

We set o''*(A) = 0. The Kondo couphng constants Ca and Cb should be in the ranges (i) Ca > 1/2, 
Cb = -3/2; (ii) Ca = 1/2, 1/2 > Cb > -3/2; (iii) 1/2 > Ca > -3/2, Cb = 1/2; (iv) Ca = -3/2, 
Cb > 1/2. For the case oi J — 2, V = —1/2, the ground state energy can be written down as the 
form 



| = ^6(0)-27ra(0), 
and (t(0) should take its largest value. Then we set p{q) — and obtain that 

'+°° b{uj) 



2M0) = 21n2+±/J- 



-exp(- \u}\) 



duj 



where 



b{uj) = S(j^^3 exp 



uj{Ca + 



Sc^-i exp 



iO Ca 



(112) 



(113) 



+Sc,+^ exp 



i^{Ca + 



S'c-i exp 



ujlCb 



exp 



(-III) 



Therefore, the finite-size correction of the ground-state energy due to impurities is 

8(2C,2 + SCa + 2) 3 TT ^^,2Ca^ 

^ = {2Ca + 3)(4C,^ - 1) - 2 - ' + 2 - 



(114) 



when Ca > 1/2 and Cb = -3/2, where /3 is defined by P{x) = 5 [?/' (^) - V' (f )] and iJj{x) = 
^ lnr(a::). By taking a ccou nt of Cg > 1/2, then Jq < 0, we have (t(A) = for the ground state. 
From relations (|lO|), ( |lO|) and ( ^T\) , we get that 



^'(A)-^ [a (a,^) +a(XCa + l)+a[X-l)-a(XCa-l) - a (x, -2) 



/(0) = 7T + 



27rC 



^-y + ln2+|-2/3(C,) 



(115) 
(116) 



p''(q) = isec/i|7rq| + ^ 



a{q,Ca + 1) -a[q, 



AttC 



+00 



6(aj) cos(ajq) 
cosh I 



for q / 0. When C^ 1/2, -3/2 < Cb < 1/2, from relation ( pISD , we have that 

'2C;, + 3 



27rcr(0) = 2 In 2 + <{ 2(1 - In 2) - TT + 2 



/3 



1 - 2C,, 



Then, the finite-size correction of the ground state has the form 
E' 



TT 








(\ 








cot 




-\^Cb) 


+ tan -TT H 




+ 1 


2 








V4 







■ln2 



1 44Cfc2 _ 26Cb - 35 40C;^ 
2(2C6 + 3)(2Cb-l)(2Ct + l)' 



(117) 



(118) 



(119) 
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By taking account of Jt > 0, we have a^{\) = 0. From relations (|l05|) , ( |l06|) and ( |l07|) , we obtain 
that 



diX) = 



1 



47rG 7o 



+°° cos(wA) 



duj- 



+Sc,+i exp 



cosh - 



(120) 



= ^sec/i|7rg| + ^ 



a [ g, - ) +a(gf,C6 + 11 



47rG 



+°° , cos(u;g) 
cosh| 



exp (— \2lu\) + exp (— |cj|) + exp 



"5, 



_i exp 



Cb 



for q 7^ and 



p\0) = 



for Cb 



l + 2;rG{c;^ + /3(^^'' + 2) + %-i/3(^^H^)} for a ^ -i 



(121) 



(122) 



The cases of 1/2 > Ca > —3/2, Cb = 1/2; Ca = —3/2, Cb > 1/2 have the similar expressions. 
When J = —2, F = 1/2, by similar discussions, the finite-size correction of the ground-state energy 
can be written down as 



E' = -- 



4(2^ + 1) 



(2a-l)(2Ca + 3) 



for Ca > 1/2, Cb = -3/2 and 



E' 



3 _ 4(2Cfc + 1) 
2 (2a-l)(2C6 + 3) 



(123) 



(124) 



for Ca = 1/2, 1/2 > Cb > -3/2. 

As the conclusions, an integrable model in one dimension is constructed from t — J model where 
two magnetic impurities are coupled to the system. It describe the behavior of the strong correlation 
electrons with Kondo problem. The spectrums of the system is not linear. The boundary R matrix 
depends on the spin and rapidity of the particle and satisfies the reflecting factorizablc condition. 
The Hamiltonian of the model is diagonalized exactly by the Bethe-Ansatz method. The integral 
equations are derived with the complex "rapidities" q which describe the bound states in the 
system. The properties of the ground state are discussed and the finite-size corrections of the 
ground-state energies are obtained due to the couplings of the magnetic impurities. 
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